Semi-microscopic description of the backbending phenomena in 

some deformed even-even nuclei 

A. A. Raduta a ' fe) and R. Budaca a ) 
a ' Institute of Physics and Nuclear Engineering, 
RO-077125 Bucharest, POB MG6, Romania and 
^Academy of Romanian Scientists, 54 Splaiul Independentei, Bucharest 050094, Romania 

(Dated: August 26, 2011) 

Abstract 

The mechanism of backbending is semi-phenomenologically investigated based on the hybridiza- 
tion of two rotational bands. These bands are defined by treating a model Hamiltonian describing 
two interacting subsystems: a set of particles moving in a deformed mean field and interacting 
among themselves through an effective pairing force and a phenomenological deformed core whose 
intrinsic ground state is an axially symmetric coherent boson state. The two components interact 
with each other by a quadrupole-quadrupole and a spin-spin interaction. The total Hamiltonian 
is considered in the space of states with good angular momentum, projected from a quadrupole 
deformed product function. The single-particle factor function defines the nature of the rotational 
bands, one corresponding to the ground band in which all particles are paired and another one 
built upon a 213/2 neutron broken pair. The formalism is applied to six deformed even-even nuclei, 
known as being good backbenders. Agreement between theory and experiment is fairly good. 

PACS numbers: 21.10.Re; 21.60.Ev; 21.60.Ev; 21.10.Hw; 23.20.js; 27.70. +q 
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I. INTRODUCTION 



The anomaly in energy spacings of the rotational spectra is still an actual subject for 
theoretical and experimental studies. A special attention is payed to the backbending phe- 
nomena observed in the moment of inertia dependency on the angular velocity squared. 
The sudden increase for the moment of inertia at intermediate and high spin is reflected in 
the energy spectra by a discontinuity of the monotonous increase in energy level spacings. 
Backbending is a common phenomenon for many heavy and quadrupole deformed nuclei. 
Since its discovery |1|], there were many attempts to provide a theoretical interpretation for 
such an anomalous behavior of nuclear energy spectra. The last decades studies showed 
that in general the backbending is a result of the crossing of the ground band (<^band) with 
another rotational band with larger moment of inertia |2H4|. The mechanism of backbend- 
ing is now well known, i.e. 213/2 two neutron quasiparticle alignments. The nature of the 
second rotational band was, however, a question of a long standing debate, such that a few 
theoretical interpretations for the origin of this band came out [5J]: i) The second band has 
a deformation which is larger than that characterizing the ground band, ii) The second 
band is not superfluid like the ground band but a rigid body one. hi) The second band is 
built upon a broken pair with aligned individual high angular momenta. The total angular 
momentum of the de-paired particles being itself aligned to the core angular momentum. 
The later two hypothesis were the most successful ones. 

The pioneering theoretical interpretation is based on the so called Coriolis Anti Pairing 
Effect (CAP) proposed by Mottelson and Valatin [6(. Indeed, the Coriolis interaction vio- 
lates the time reversal symmetry and consequently contributes essentially to the de-pairing 
process. As a result, a phase transition from a superfluid to an independent particle state 
characterized by normal fluid properties and therefore by a larger moment of inertia, takes 
place. CAP is analogous to the Meissner effect in metal superconductors [7]. The Coriolis 
force is proportional to the orbital angular momentum I, which results in breaking first the 
pairs built on nucleons having largest I with respect to the axis of rotation. Breaking a pair 
leads to a dramatic increase of the moment of inertia [l| due to the large decrease of the 
static gap which may even vanish. Within a cranked mean field approach, the backbending 
phenomenon is caused by a rearrangement of the vacuum configuration or alternatively by 
the crossing of the ground-state band with the lowest two quasiparticle (2qp) band which 
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is often referred to as the S'(tockholm) band jlj. Information about several features like 
the crossing frequency, the yrast-yrare interaction and the alignment gain can be obtained 
from the diagram showing the qp Routhians calculated at fixed deformation and for a con- 
stant pairing gap, versus the frequency 8KL3J. Difficulties raised due to the semiclassical 
nature of the cranking approach and related with the angular momentum dispersion, A(J 2 ) 
were discussed by several authors (see for example 12]). Certainly a quantitative micro- 
scopic description of the moment of inertia in the crossing region encounters some difficulties 
caused by the symmetry breaking. Thus the Lipkin-Nogami method [3, [ijj for the gauge 



)rojection was used in Ref. 



15] by Satula and collaborators, while the Galilean invariance 
16] has been restored in Refs.[l3, Q- The spurious shape dependence met in treating 



the q uadr up ole- q uadrupole (QQ) pairing interaction has been eliminated by using a double 



stretched QQ pairing 



15] and an excellent quantitative description of the moment of iner- 



tia for the so called super- deformed bands of Hg-Pb nuclei was obtained. An alternative 
description of the pair breaking process [19] was proposed by Stephens and Simon, pointing 
to a mechanism which causes a rotational alignment of the particles from intruder orbitals. 



Other microscopic models reproducing qualitatively the zigzags 



lape of experimental plots 



of moment of inertia are those of Sorensen 20] and Faessler 21, 22]. More recent attempts 
proposed approaches based on the Interacting Boson Model 23l. |24|. 

In this study we present a new and simple semi-phenomenological model for backbending 
which here is considered to be the result of the band crossing mechanism between the ground 
band and a two quasiparticle decoupled band known as the 5-band. The first backbending 
is known to be caused by the breaking of a neutron pair from the intruder orbital £13/2, 
while the second one is due to a subsequent breaking of a proton h n / 2 pair. Our model 
is meant to reproduce the low spin and intermediate spin states from the yrast band, such 
that the scope of this paper is to describe only the first backbending. The particles from the 
intruder orbitals where the pair breaking occurs are treated separately from the remaining 
ones which define the core. The phenomenological core is described by the Coherent State 
Model 25] , while the motion of the intruder particles is treated through BCS model states. 
A special ingredient of our formalism is that we consider here a deformed core which induces 
also deformed trajectories for the intruder particles. The rotational bands implied in the 
hybridization procedure are defined by angular momentum projection from quadrupole de- 
formed product states, which achieves the coupling of the single-particle degrees of freedom 
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to a quadrupole deformed core and provides states with good angular momentum. The basis 
states are deformed and non-orthogonal. From this basis an orthogonal basis is obtained 
which is further used to diagonalize the model Hamiltonian. The yrast energies are defined 
by the lowest eigenvalues of the model Hamiltonian in the orthogonal basis. This formalism 
was applied to six even-even nuclei from the rare earth region with N = 90 — 94 neutrons, 
which exhibit backbending behavior in their moment of inertia plots. The numerical appli- 
cations reproduce quite well the experimental data and also provide some useful information 
regarding the rotational alignment of particles moving in the intruder orbitals. 

Results are presented according to the following plan. In the next chapter we describe the 
projected non-orthogonal basis states from which emerges the orthogonal basis. The model 
Hamiltonian is introduced in Sec. Ill where some of its single-particle and collective features 
are discussed. The matrix elements of the model Hamiltonian in the non-orthogonal basis 
are analytically given. Sec. IV is devoted to deriving the final orthogonal basis and to the 
band hybridization procedure. Numerical results are analyzed in the Section V. Conclusions 
and the future perspectives of the model are discussed in Section VI. 

II. DESCRIPTION OF THE BACKBENDING PHENOMENA IN A RE- 
STRICTED MODEL SPACE 

As we have already mentioned a particle-core interacting Hamiltonian will be treated in 
a space of spherical projected particle-core states. We start by giving the necessary details 
for the projection procedure and therefore about the model states. The collective factor 
state is a coherent state for the quadrupole bosons b\ , while the single-particle component 
is a deformed BCS state describing a set of paired nucleons which are moving in a deformed 



mean field 



26|]: 



= ^ 6 = \BCS) d e d ^- b ^\0} b . (2.1) 

In the above equation, |0)& stands for the boson vacuum state, and d is a real parameter 
which simulates the nuclear deformation. The low index of the BCS state suggests that this 
is deformed. The projected states are obtained, in the usual manner, by acting on the state 
(2.1) with the Hill- Wheeler projection operator 



P mk = / Dti K R(Sl)d£l. (2.2) 
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In what follows, we shall present the angular momentum projection procedure for each factor 
state from (2.1). 



A. Angular momentum projection of the coherent state 

The projection of the quadrupole coherent state was presented by one of authors (A. A. 
Raduta), in collaboration, in Ref. 27J. The main results are as follows: 

= UfP J M ^ b . (2.3) 

The matrix elements of any boson Hamiltonian between these projected states can be ana- 
lytically expressed in terms of the norms: 

(A4 s) )~ 2 = (2J + l)e~ d2 lf\ (2.4) 

where Ij stands for the overlap integrals 

if = [ Pj(y) [P 2 (y)] k e* p ^dy, with x = d\ (2.5) 
Jo 

Pj{y) denotes the Legendre polynomial of the rank J. These integrals have been analytically 
calculated in Ref.Q, 0. 

B. Spherical projected states from a deformed BCS state 

The projection procedure adopted in this paper for the deformed BCS state is that 
formulated by Kelemen and Dreizler in Ref. |28| . Thus, the fermionic function ipf is expressed 
first as a linear combination of states with a definite angular momentum 

1> f = \BCS) d = Y,Cj\J,0), (2.6) 
j 

and then the projection operator selects only the component of the desired angular momen- 
tum: 

= M { j f) P^ f = Afj f) Cj\J, M). (2.7) 
The projected function is normalized to unity, and consequently we have: 

Afy ) = Cj\ (2.8) 
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In this way the amplitude Cj can be expressed as: 

2J + 1 



\ C J2_ U. loJ 



&s\PMf) 



P J (cos/3)(^ / |e- i/3 ^|V/) sin 0dP 



(2.9) 



The calculation of Cj is very much simplified if the projector operator Pq is expressed as a 
finite sum of particular rotation operators [281 ]: 



oo 



2J+1 
M + l 



M/2 



71=1 



(2.10) 



where J is even and M = J max - On the other hand the terms Ai(M, J) and B n (M, J) can 
be also analytically calculated. Their expressions are given in Appendix A. 

The calculation of the amplitudes Cj is thus reduced to finding the matrix element 



^ nJv \ipf). For this purpose, the fermionic wave function is written as a sum of 



components with determined number of particle pairs: 



\Bcs) d = J] (u 3m + v 3mC ) m c)_ m {-y- 

m>0 

,m>0 / \ fci>0 ]Kl 



10) 



E 



fci<fc 2 



r/-, 77., * ; u jk 2 u j-k 2 \ J 



+ . . 



(2.11) 



where \N P ) are states with N p pairs of particles and {U, V} are defining the Bogoliubov- 
Valatin (BV) transformation from the particle to the quasiparticle representation. It can 
be shown that the matrix element of the rotation operator e~ lJT + I ' n " Jy on states \N p ) can be 
expressed in the form of a determinant of rank 2N P of reduced Wigner functions d J ki k2 (/3^ p ) 
with the argument: 

"""> = WTi' (2 ' 12) 

where M is the maximum angular momentum realized by the set of N p pairs of particles. 
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For illustration we will present here only the case for two pairs: 



(o|(- 



c j—ka c jka \ 



I ^j-ki^jkiZ jk[ j-k 



' _){ki+k 2 +k' 1 +k' 2 ) ^ 



(2.13) 



i(-)^ fci i^-^(-)^ fc Ho) 

/ di iK m di u _ k ,m 4 uk ,(^) d? ku _ k ,m ^ 
d\ K m di ku - k[ m JL^m di ku _ k ,m 
<, k[ m di,- k[ m di k ,m di 2 _ k ,(^) 
\di k2K m) di k2 ^,m) d>_ k2A m d^_ k ,m j 

The largest angular momentum M in the configuration (j) 2Np , where j is the angular mo- 
mentum of the individual particles, is not necessary 2N p j because of the Pauli principle 
constraint. Group theory provides a simple formula for the upper limit of the total angular 
momentum of a given configuration, which moreover takes care of the Pauli principle [29] : 



M = N p (2j - 2N P + 1) 



(2.14) 



In Table I the largest angular momenta achieved for all possible numbers of neutron pairs 
in the intruder orbital ^13/2, are listed. The results from this table are used to calculate the 



TABLE I: Maximum angular momentum M achieved for a given number of particles, each carrying 
an angular momentum of j = 13/2. 



Number of particles, 2N P 





2 


4 


6 8 


10 


12 


14 


Maximum angular momentum M 





12 


20 


24 24 


20 


12 






arguments by means of Eq. fl2.12p and then the overlaps of the type (I2.13p . With all 
these done the average of the rotation operator e~ l M+\ nJ y with \BCS)d is readily obtained. 
Finally the projected particle-core function is written in the form: 



(i) 

JM 



jWPito\BCS)rt a =N?>Y, c ^ 



, j Nf cs N { j 9) 

JfJc Jj Jc 



J JM 



with the normalization factor 



a 



JfJcJ 





, , \Nf cs N[ 9) 

J f Jc \ *J f Jc 



(2.15) 



(2.16) 



The summations in Eqs. (I2.15P and (I2.16P are restricted to the ranges Jf < 24 (see Table I), 
and J c < 60. In this way one accounts for all possible configurations of a given total angular 



momentum J, with J running from to 36. The set of wave functions (12 . 1 5[) describes the 
ground band which is associated to the case when all particles from the intruder orbital are 
paired. 

The 5-band which crosses the <^band and produces the backbending in the energy spectra 
can be described by a wave function similar to ( I2.15P with the difference that now one pair 
of particles is broken, i.e. they occupy two states which are not related by a time reversal 
transformation. The symmetry breaking is simulated by applying the angular momentum 
raising operator on a function with good symmetry. Thus the two qp state which is respon- 
sible of generating the 5-band is a K = 1 state of the form: J + aj fc aj_J.BCS%, where aL 
is the creation quasiparticle operator defined by the canonical BV transformation: 



U 3kc) k - Vj 



k[, )"' Cj— An 



) j - k ct 



OCjk — vjk^jk y jk\ ) ^j-k- 

The total projected state corresponding to the S^band is defined by 



(2.17) 



J 

MX 



J + a) k ay k \BCS) c 



JfJc 



CJ f J C J 
1 1 



jk As) 



JM 

(2.18) 



where the normalization factor is 



-2 



E 

JfJc 



CJ f J c.J 
1 1 

A^; i .v.r 



(2.19) 



The projected two quasiparticle state &j \ from (I2.18P has the following expression: 



jk 



■J/1 
P k TD J f 



■j f i;M f ~ AX: :•/•'>,,'>; , "C*h- 
The normalization factor of this state is defined by the matrix element 



Af- 



■jk 



-2 



i(BCS\aij k aij-kJ-Pii J+oi] k a\_ k \BCS) d , 



(2.20) 



(2.21) 



which is determined in a similar way as the norm of the projected BCS function (see Ap- 
pendix B). The upper limit of the Jf is still 24 even after the application of the quasiparticle 
creation operators, because it is the maximum angular momentum which can be realized in 
the i\3/2 intruder orbital, irrespective of the number of particles involved. 



Both projected states, (I2.15P and (I2.18p . depend on the deformation parameter d, al- 
though they are states with good angular momentum. Note that the projected 2qp states 
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are defined only for even angular momentum J, with J > 1. Another important property 
of the projected 2qp states consists of that two states with k ^ k' are not orthogonal. 



III. THE MODEL HAMILTONIAN 

The backbending features of some rare earth nuclei will be studied with the following 
particle-core Hamiltonian: 

H = H c + H, p + H pair + H pc . (3.1) 
The Hamiltonian H c is a harmonic quadrupole boson operator: 

H c = hu b J2bl^, (3.2) 

and describes a spherical core. H p is describing a set of particles in a spherical shell model 
single j orbital of intruder nature, interacting through a paring force: 

Hp = (s n lj — A) C nljm C nljm, 
m=all 

G 

T 



where Pj(Pj) are creation(annihilation) operators of the Cooper pair in the intruder orbital 
j defined by 

p} = T, c Wr.(-r ra ( 3 - 4 ) 

rre>0 

The operators - m and c n ^ m stand for the creation and annihilation operators for a particle 
in the spherical shell model state \nljm) having the energy e n ij. The Lagrange multiplier A 
plays the role of the Fermi energy for the paired system. For the sake of saving the space in 
what follows the spherical shell model state will be specified only by two quantum numbers, 
i.e. \jm). 

The particle-core interaction is taken of the form: 
Hp C = H qQ + H JfJc 

= -^C J2 mm ') C )m C jm' ((-r4- M + ^2 M ) + CJ C ■ Jf, 

q 2fJ ,(j;mm') = (jm\r 2 Y 2ll \jm') , (3.5) 
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where Ac and C are free parameters. The last term plays an important role in reproducing 
the correct transition at the band crossing point and simulates the effect of the Coriolis 



coupling [30 ]. 



Depending on whether the considered system is a near spherical or a deformed nucleus we 
diagonalize first the single-particle plus the pairing interaction and then treat the remaining 
terms or we diagonalize first the single-particle plus the qQ interaction and then treat the 
rest in the resulting single-particle basis 3l|. Actually, here we make the option for the 
second procedure for reasons which will become clear in the next subsection. 

A. Pairing in a deformed single-particle basis 

The particle-core interaction leads to deforming the single-particle mean field. On the 
other hand the interaction deforms the quadrupole boson. This mutual deformation effect 
is suggested by the following reasoning. If one considers the average of the particle-core 
Hamiltonian with a single-particle state one obtains a deformed Hamiltonian with a ground 
state described by an axially symmetric coherent state. As for the single-particle deformed 
mean field let us consider, for simplification, the model Hamiltonian which corresponds to 
vanishing spin-spin interaction i.e., C = and unpaired particles (G = 0), 

H = H p + H c + H q q. (3.6) 

Averaging this Hamiltonian with the coherent state 



V>& = exp 



d(bl -b 20 )]\0), (3.7) 



one obtains a single-particle Hamiltonian for a deformed mean field similar to that used 



within the Nilsson model 



22| : 



Hp = d 2 ftw b + Hp- 2dA c ^ <?2o(j; m, -m)c) m c jm . (3.8) 

rn 

Note the role of nuclear deformation played by the parameter d. Apart from an additive 
constant, in the first order of perturbation, the energies of H p are given by: 

enijm = e nlj - AdX c (2n + 3)C{ \ {C^i, (3.9) 

2^2 
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where 

Here M and Uq are the nucleon mass and the frequency of the harmonic oscillator function. 
The dependence of single-particle energy on deformation is linear and shown in Fig.l for the 
n = 5 shell. 

We remark that in a single j calculation, the n and I quantum numbers are superflue 
and therefore dropped out. Moreover, since only relative energies are involved in the BCS 
calculations the constant term, i.e. that one not depending on m, is put equal to zero. 
Therefore, in our calculations the single-particle space is spanned by \jm) and the single- 
particle energies corresponding to the mentioned states are: 

e jm = -4dX c (2n + 3)C{ 2 JCili, (3.11) 

2^2 

with n being the principal quantum number characterizing the major shell to which the 
intruder state belongs. The pairing interaction in such a deformed multiplet has been 
considered by Bes at al in Ref. [26 1. Note that the deformation effect due to the qQ 
interaction onto the single-particle state \jm), was ignored at this stage. According to Eq. 
(13.1 ip the energies of the time reversed states are equal. Therefore, we can restrict the space 
to the states \jm) with m > keeping in mind that on each such state two nucleons are 
allowed. 

The single-particle Hamiltonian for the deformed mean field can be written as: 

= £ jmc] m C jm- (3-12) 

m=all 

The next step is to treat + H pa i r through the BCS formalism, using the BV trans- 
formation (12.171) . In the quasiparticle representation the mentioned Hamiltonian becomes: 

H qp = E + ^ E ]m a )m a im + ^ 5 <jU~} ^" '("jm"] -m + "j-mC^m), (3-13) 
m=all m>0 

where the following notations were adopted: 

F - \N A ' V ~ A (^n-A) + A 2 e jm A 

i^o - -AiV part - — , ti jm - , g jm — , [Ej - U). (a. 14) 

Here N part represents the number of particles in the intruder orbital, Ej m - quasiparticle 
energy while A is the energy gap. 
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SfS/2 5/2 6i 13/2 13/2 



5p 3/2 3/2 



3 



to- 




PIC 1: (Color online) Neutron energies of deformed single-particle states given by Eq. (3.9) 
with Xc = 50.81 keV are plotted as function of the deformation parameter d. The solid thick 
lines correspond to the i 13 / 2 states. Energies are given in units of hui (= AlA'^MeV). The 
intersection of the dashed line, d = 3.103, with the energy curves indicates the single-particle 
energies which correspond to 160 Er. 

The expression of the qQ interaction term in the qp representation is: 



qQ 



2 mm ) V jm + Q20U; mm)(Uj m - Vf m )a] 



m>0 m=all 



4o + fo 20 . (3.15) 



-2 Y <?2o(j; mm)U jm V jm (-Y m (a) m a) _ m + aj_ m a jm 

m>0 

In deriving this expression we took into account the fact that only the component with 
/i = of the boson factor contributes when the average on the projected coherent state is 
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calculated. 



B. The diagonalization of H in the particle-core product basis 

We recall the fact that the BCS equations can be obtained either by minimizing the 
ground state energy or by vanishing the dangerous graphs. Indeed one can check that the 
total coefficient multiplying the cross terms (— y~ m (oij m aij m + aij - m aj m ), coming from H qp 
and H q Q is vanishing. 

The diagonal matrix elements of the various terms of the model Hamiltonian can be easily 
calculated. Their compact expressions are: 

(¥p M \H qp + H gQ \¥p M ) = E + 2dA eff , (3.16) 



(*f M Jjk)\H qp + H qQ \^? (jk)) = E + 2E jk + 2dA eff , (3.17) 



where 



A eff = -2X c (2n + 3)C{H £ C^Mr ( 3 - 18 ) 



2 " 2 

m>0 



is a constant quantity and can be omitted along with Eq. 

The matrix elements of the harmonic boson Hamiltonian (3.2) on projected coherent 



states were given in Ref. 33], where a description for the ground band energies of axially 
deformed nuclei was provided. Thus, the matrix elements of H c in the space of Oqp and 2qp 
states are given by: 

JfJc \ Jf Jc / Jc 



2 



{^■m\ H A,m) = « 2 (Kw) 2 Yl (■0jfi) $)■ (3 - 20) 

JfJc \ J f^- J c f J c 



Using the tensorial form of Oqp and 2qp states (12.151) . (I2.18P as well as the results of 
Appendix B, the matrix elements of the spin-spin interaction term Hj f j c are easily found: 

<*ai^i*a> - ^ + i)^K i ) ! E( s i|y) 1 

JfJc \ Jf Jc / 

x [J c (J e + l) + Jf(J f + l)], (3.21) 

^%-m\HjfjAm) = (M?w) 2 £f J^L) 

JfJc \^Jfl N Jc J 

x [J c (J c + l) + J f (J f + l)]. (3.22) 
13 



In this way the rotational spectra of g^band and 5-band are completely determined. 



IV. BAND HYBRIDIZATION PROCEDURE 

To describe the backbending phenomenon as being determined by the bands intersec- 
tion we have to diagonalize the total Hamiltonian in a basis defined by (]2.15j) and ( 12. 18ft . 
Unfortunately, this basis is not orthogonal. Indeed the overlap matrix elements 

0( 2 M (jk) = (*?iO'*)l*Sf;lW>, (4-1) 

are not vanishing. This can be seen from the explicit expression: 

J f J c iV Jc JV J/ ■ /V J/1 

Overlaps are nonvanishing due to the fact that quasiparticle operators are not tensors of a 
definite rank and projection, which can be easily seen from their transformation against an 
arbitrary rotation R: 

X^kR- 1 = £ Dl k a) m + ^(-r m «i-m, (4-3) 
m 

D Lk = ( U jkUjm + VjkYjm)D 3 mk , 
D Lk = (UjkV jm - VjkUj m )D J mk . 

That would not happen if the U and V coefficients were not dependent on the m quantum 
number. By diagonalizing the overlap matrix, one obtains the eigenvalues a m (J,jk) and the 
corresponding eigenvectors Vn m \j,jk). Then the functions 

= [a m (JJk)r 1/2 [*%vl m) (J,jk) + *f M . tl Vt m \jlJk)] , m = 1,2, (4.4) 

are mutually orthogonal and can be used as a diagonalization basis for the total Hamiltonian 



34|. If we take the total wave function of the form 

m 

then we have to solve the following eigenvalue problem 

H mm ,X^(jk) = E JM (jk)Xi M (jk), (4.6) 

ml 
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for finding the yrast spectrum. The Hamiltonian matrix H nm is defined by 

#w = [«_(^ jA:)^^^, jA:)]- 1/2 ^ ^A:) <^ l^l^? )^ jA:), (4.7) 

nn' 

where the low indices /?,/? are either JM if the corresponding upper index n (or n') is 1 or 
JM; 1 for n (or n ) equal to 2. By mixing the non-orthogonal states in the orthogonalization 
process and then diagonalizing the model Hamiltonian a natural interaction of the primary 
bands takes place. Such an interaction is angular momentum dependent and seems to be 
more efficient in the K = 1 band hybridization process 35] than a constant one. 

The off-diagonal matrix elements of the model Hamiltonian terms, in the non-orthogonal 
basis are given by the following expressions: 

(*%\H qP +H qQ \¥? M;l m = 

x {E + E jk + 2dA eff ) , (4.8) 



2 U^n™**, 

x (J(J+l)_ j c (j c + i)_ J^Jz + l)). (4.10) 



V. NUMERICAL APPLICATION 

The formalism described in the previous sections was applied to six even-even nuclei 
from the rare earth region which are known to be good backbenders, namely 156 Dy, 160 Yb, 
158 ' 160 Er and 164 > 166 Hf. First three nuclei are N = 90 isotons, next two are N = 92 isotons 
and the last one has 94 neutrons. 

The model Hamiltonian involves four free parameters namely, the pairing constant G, 
quadrupole-quadrupole interaction strength X c , the spin-spin interaction strength C and 
the boson frequency of the core hwt,. Another parameter is the deformation parameter d, 
defining the coherent state 

The input data for the BCS equations are the pairing constant G and the single-particle 
energies determined by Eq.(3.9). The deformation parameter d and the boson frequency fkJb 
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TABLE II: The pairing strength (G) is given in units of MeV, while the quadrupole-quadrupole 
(Xq) and spin-spin (C) interactions strength are given in units of keV. The boson frequency of 
the core is also given. The list of the deformation parameters d are presented together with the 
corresponding /?2 deformation, taken from Ref. 



36]. The manner in which these parameters were 



fixed is explained in the text. 



Nucleus 


h 


d 


hw b {MeV] 


X c [keV] 


G[MeV] 


C[keV] 


156 Dy 


0.211 


3.1414 


0.9781 


52.63 


0.1786 


72.27 


i60 Yb 


0.195 


2.5000 


1.0078 


65.59 


0.1854 


36.90 


158g r 


0.203 


2.7355 


0.9778 


59.95 


0.1814 


60.10 


160g r 


0.231 


3.1030 


1.0089 


50.81 


0.1769 


93.63 


164 Hf 


0.208 


2.6605 


1.1039 


59.89 


0.1836 


52.20 


166 Hf 


0.237 


2.8505 


1.0436 


52.89 


0.1737 


69.63 



are fixed so that the first energy levels lying before the band crossing point are reproduced 
with a good accuracy by Eq.(3.19). Indeed, the g-band energies are not very sensitive to 
the single-particle degrees of freedom. Practically, the only contribution to the total energy 
is due to the core, the particle-core interaction being reflected in the norm of the projected 
BCS state. In the case of the g-band states, all particles are paired, and the dominant term 
of the sums in Eqs.(3.19), (3.21) corresponds to the situation Jf = and J = J c . This is 
also consistent to the fact that until the band crossing point the whole angular momentum 
is carried by the core. In this way the g-band energies can be roughly approximated by: 



E, 



]2_J 
I 



hw b d 



(5.1) 



./ 



which isiust the expression of the ground band energies predicted by the coherent state 
model {25I, 0]. In the extreme limits of large and small deformations d, the above energies 
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have been approximated by compact and simple functions of J(J + 1) in Refs 
Actually, this simple expression can be used to determine the parameters d and froj b by 
fitting the first energies before the band crossing of the g and S bands. Later on a tuning fit 
can be achieved by using the eigenvalues of the model Hamiltonian in the orthogonal basis. 

The pairing interaction constant G and the qQ interaction strength Xq are fixed such that 
the band crossing point and the observed sequence of single-particle energies are reproduced. 
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Here we deal with neutrons from the shell 113/2, where we can expect at most 7 pairs. The 
number of neutrons considered out of the core, in the shell 213/2, is 2 for 156 Dy, 158 Er, 160 Yb, 
164 Hf, 160 Er and 4 for 166 Hf. Solving the BCS equations one obtains the gap parameters 
A, the Fermi level energy A and consequently the occupation probability parameters U 
and V. With this data the wave functions corresponding to Oqp ( 12.1 5p and 2qp ( 12.1 8p 
are completely determined. The relevant information yielded by the BCS calculations are 
presented in Table III. We mention the fact that our choice for the number of neutrons 
distributed on the i 13/ / 2 sub-states is consistent with the BCS calculations in an extended 
single-particle space. Indeed, we solved the BCS equations for a space of 23 states, consisting 
in the union of the sub-states 5/i 9 / 2m , 5/ 7/2m , 6ii 3/2m , 5p 3/2m , 5/ 5/2m , 5p 1/2 m, 5h n /2,n/2, where 
we distributed 10,12 and 14 neutrons for the N = 90, 92, 94 isotons respectively. The results 
for the fitted pairing strengths can be interpolated by the following function: 

G = j (g - - with go = 37.44MeV, gl = 62.25MeV. (5.2) 

Summing up the occupation probabilities for the zi 3 / 2 sub-states, we obtained the average 
number (iV**^ 2 ) listed in Table 3. Thus N part is the largest even number smaller than (N^t)- 
Exception is for the case of 166 Hf where N part is larger than, but very close to, (iV „% 2 ). The 
calculations in the reduced space of ii 3 / 2 sub-states was performed with a pairing strength 
chosen such that the minimal qp energy be equal to that obtained within the extended 
space. In this way the occupation probabilities obtained by solving the BCS calculations 
in the extended and reduced single-particle space respectively, are close to each other. In 
his context one could say that the schematic calculations in the reduced single-particle space 
accounts for the effective pairing interaction in this subspace. 

Note that except for 166 Hf, the gap A is very close in magnitude to the corresponding 
quasiparticle energies, which suggests that the Fermi level lies close to the selected /c-energy 
level. Of course, the S"-band associated to the 2qp projected state (fc) 2 is the first one which 
intersects the (/-band and will become the yrast after intersection. The quantum numbers k 
listed in Table III are consistent with the Nilsson model prediction for the last filled orbital 
of the chosen nuclei having the quadrupole nuclear deformation /3 2 given in Table II. 

Our calculations show that the Jf = component has a squared weight of about 45-64% 
in the BCS composition. As a matter of fact, this suggests that the approximation (5.1) 
made for the <^band energies is to be corrected due to the nonvanishing angular momentum 
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TABLE III: The Fermi level energies, gap parameters and the quasiparticle energies are given for 
the chosen number of particles N par t, and the projection k associated to the broken pair. We also 
give the average number of particle, denoted by (iV*^ 2 ), in the multiplet ^13/2- 





N 

1 V 


part 


(N h3/2 \ 

\ part 1 


a 1 ivie v 1 


A (MpVI 
i-i 1 ivie v 1 


lei H,i TT1 1 tl I 

fb I J-Jfc 1 1 bit L \ 


F,i nVTpVl 


156 Dy 


90 


2 


2.39 


48.7643 


1.12446 


1/2 


1.12448 


160 Yb 


90 


2 


2.47 


48.3187 


1.19819 


1/2 


1.19830 


158 Er 


90 


2 


2.44 


48.5354 


1.15443 


1/2 


1.15460 


160 Er 


92 


2 


3.12 


48.6732 


1.16028 


3/2 


1.16036 


164 ffi 


92 


2 


3.18 


48.2444 


1.23826 


3/2 


1.23827 


166 jj f 


94 


4 


3.88 


48.4181 


1.17327 


3/2 


1.20183 



components of the BCS state. In fact since the Jf = component is only moderately 
dominant in the BCS state composition, the effective angular momentum of fermions is not 
vanishing in the g band and not 12 for the states belonging to the S band. This feature is 
conspicuous from the analysis of Fig. 5. 

The matrix elements between states (2.18) of the spin-spin interaction term have a very 
peculiar feature. Their numerical values increase with the total angular momentum, going 
from negative values to positive ones. The transition from negative to positive values takes 
place at J = 10, 12, where, in most cases, the backbending shows up. In this situation the 
parameter C does not affect the position of the band crossing point. However, it influences 
the magnitude of the yrast energies in the high spin region. Due to this feature C is fixed 
as to reproduce the moderate high spin energies of the yrast band. The fitted parameters 
of the model Hamiltonian are collected in Table II. We remark that the factor dXc which, 
according to Eq.(3.9), plays the role of the deformed mean field strength, has roughly a 
linear dependence on the quadrupole nuclear deformation /3 2 - This dependence is illustrated 
in Fig.2. 

Using the parameters specified in Table II the equation ( 14. 6 p provides the system energies. 
For each angular momentum J we select the lowest eigenvalues E(J). The set of energies 
Ej defines the yrast band. The discrete derivatives of the resulting energies with respect to 
the angular momentum defines the angular frequency: 

M J) = ^dP- W I [E{J + 2) ~ E{J)] ■ (5 ' 3) 
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FIG. 2: (Color online) The product dXc, with the factors given in Table II, is represented as a 



function of the nuclear deformation 02 taken from Ref. 



Alternatively, the angular velocity can be also denned by using for E(J) the expression 
provided by a symmetric rotor Hamiltonian: 

E(J) = (5.4) 

Then the discrete derivative of this expression yields: 

MJ) = ^j*- (5-5) 

From here one derives a simple expression for the moment of inertia: 

4J + 6 

J ~ E(J + 2)-E(JY [b - b) 

The backbending plot is a graph in which the moment of inertia is plotted versus (ftw) 2 . 
Theoretical results and experimental data are usually compared in terms of this plot. This 
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is done in the Fig. 3 for the even-even rare earth nuclei treated here. Note that in all 
cases the zigzag behavior is reproduced quite well. In general, a good agreement between 
theoretical and experimental data corresponding to low-spin states and moderately high 
spin states located however below the possible second backbending is obtained. The second 
backbending is known to be caused by the consecutive breaking of a neutron i 13 / 2 pair and 
a proton pair from the shell hn/ 2 . The second backbending seems to be a rare event. 
Despite of this some of the considered nuclei exhibit such a phenomenon. Indeed, 160 Yb has 
a second backbending which occurs at J = 26, while for 158 Er the second backbending shows 
up in the same region but is less evident being rather an up-bending. In both mentioned 
cases our calculations show an up-bending around the J=26 state. The study of the second 
backbending is however beyond the scope of the present paper. 

The crossing of the g and the S bands is illustrated in Fig.4, where their energies are 
plotted as function of J( J+l). We also give the results for the two unperturbed bands whose 
energies are approximated as the diagonal matrix elements of the system's Hamiltonian in 
the non-orthogonal basis. In each panel the r.m.s. values for the deviation of the theoretical 
results from the corresponding experimental data are given. As can be seen from Fig.4, 
while the 5-band energies exhibit a linear dependence on J(J + 1) for the g-band a slight 
quadratic dependence on J (J+l) is met. Since 156 Dy and 160 Er have the largest deformation 
d, the linear dependence on J(J + 1) of the corresponding g-band energies prevail. 

In order to investigate the alignment of the individual angular momenta to the angular 
momentum of the core we define the average angular momenta for the interacting compo- 
nents as follows: 

1(1 + 1) = (* J £m\E\* J £W), (5-7) 

j f (j f +i) = (® j T %m\jJ\®Tot(m- (5.8) 

The full alignment corresponds to the situation when J c + Jf equates the total angular 
momentum of the system J. The departure from this ideal picture is measured by the 
deviation: 

AJ = J-(J C + J f ) . (5.9) 

This quantity together with the average angular momenta of the core J c and fermion system 
Jf are plotted versus the total angular momentum J, in Fig. 5. These plots reveal several 
features, such as the band crossing point, the amount of angular momentum carried by the 
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FIG. 3: (Color online) Backbending plots for 156 Dy, 160 Yb, 158 > 160 Er and 164 > 166 Hf isotopes com- 
paring theory (black squares) with experiment (red circles). Experimental data are taken from 
3J-|43]. 
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FIG. 4: Energy trajectories resulting from the diagonalization of total Hamiltonian in the orthog- 
onal basis (4.4) (straight lines) are compared with experimental points (black circles) for 156 Dy, 
160 Yb, 158 ' 160 Er and 164 > 166 Hf. The unperturbed #-band (dotted lines) and the selected 2qp band 
(dashed lines) are also presented. For each nucleus the r.m.s. values corresponding only to the 
states taken in consideration here, i.e. up to J = 26, 28, are given. 
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broken pair or even the fraction of the angular momenta alignment. For example for Er 
isotopes and 166 Hf the band crossing takes place at J = 12, while for 164 Hf and 160 Yb at 
J = 10 and for 156 Dy at J = 14. This is actually a confirmation of the backbending results 
from Fig. 3. The amount of angular momentum, Jf, carried by the broken pair varies from 
10h to lAh. We remark that at the band crossing the alignment defect is maximal and this 
is decreasing by increasing J. The full alignment is never reached but the defect A J exhibits 
a plateau with the value of 2-3h, beyond J = 20. This is a salient feature for the present 
formalism which is not met in other approaches where AJ C = A Jf = 10, 12. Indeed, for 
160 Yb the variation of J c around crossing is only 2h and, moreover, after band crossing the 
core and fermionic angular momenta are almost equal to each other. For all other nuclei, 
at the crossing J the core angular momentum variation is about Ah. We don't have any 
situation where after crossing J c = . Therefore for a critical value of angular momentum 
of the core the spin-spin interaction causes the de-pairing of two neutrons, which almost 
align their angular momenta and starting from a larger total spin (~ 20) the fermion total 
angular momentum is almost aligned to the core angular momentum, AJ = 2 — 3h. In 
general, the second alignment is produced for J c > Jf 

Figure 6 shows the calculated rotational energies and the experimental data of g-band 
and 5-band versus angular momentum for 166 Hf. We present this case separately since here 
the best agreement of theoretical results and the corresponding data was obtained in the 
backbending region, after the crossing point. The slope of the curves from this figure is the 
angular velocity ou(J). Note that up to a critical spin the rotational energies of the 2qp band 
have a negative slope then this is vanishing and further is increasing. Finally it reaches 
a constant value which reflects an equidistant structure for the spectrum in this region. 
The later situation shows up when a maximal alignment is achieved. The negative slope 
corresponds to the situation when the fermion and the core angular momenta make an angle 
varying from n to ir/2. At the beginning of the interval the individual angular momenta are 
almost anti-aligned and their sum is also almost anti-aligned to the core angular momenta. 
When the individual angular momenta are partially anti-aligned, the two frequencies caused 
by particles and core respectively, add destructively. Increasing J the two frequencies have 
the same sign and therefore add each other, constructively. In this part of the curve, the 
second alignment namely that of the total fermionic and the core angular momenta, starts 
operating. When a maximal alignment is achieved the slope keeps constant and consequently 
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FIG. 5: (Color online) Expected angular momentum of the core and of the broken pair of intruder 
neutrons. The deviation of total angular momentum A J is also presented. 

the two curves become almost parallel. As shown in Fig. 5, the properties described before, 
seem to be generally valid. A similar angular momentum dependency was obtained also in 
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FIG. 6: (Color online) Energy plot for 166 Hf showing the interacting bands. Rotational energy of 
the g-band and of the 2qp band from a certain spin increases monotonously. Experimental data 
points are also visualized. 
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44|. 



The mechanism of breaking the pairs and aligning the individual angular momenta to 
the core angular momentum is suggested in Fig. 7. Indeed, for low values of the core angular 
momentum the nucleon angular momenta are aligned to the symmetry axis of the mean field 
which, as a matter of fact, is determined by the qQ coupling term. This situation is shown 
in Fig. 7 a) where J/ = and J c = J. 

The mechanism suggested in Fig. 7 is consistent with the following phenomenological pic- 
ture. The core angular momentum is perpendicular to the symmetry axis OZ. Let us consider 
that OX is the direction of J c . Averaging Hj f j c with the Wigner functions D J ^ K with M 
and K being eigenvalues of (J c ) x , one obtains an operator acting in the intrinsic frame which 
breaks the time reversal symmetry. Due to this property this is the term responsible for 
the neutron pair breaking connecting the states of Oqp and 2qp types, respectively. In the 
intrinsic frame the mentioned term connects a state with all particle paired, i.e. with total K 
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equal to zero, with a state where a K = pair is replaced by a K = 1 broken pair. Actually, 
the same effect is obtained within microscopic models using for a many body system with 
pairing interaction a cranking term of the type uj x . It is worth mentioning the following 
specific feature of the present model. The total quantum number K is obtained by summing 
the contributions from nucleons and the core. However, writing the projected state of the 
core [25] in the intrinsic frame, one obtains an expression which is a superposition of com- 
ponents of different K. Among these, the major component corresponds to K = 0. In this 
respect the broken pair state is not a pure K = 1 state but, due to the core, a superposition 
of various K ^ components. However the major component is the one characterized by 
K = 1. 

At this stage it is worth making clear the way the spin-spin interaction simulates the 
action of the Coriolis force in the intrinsic reference frame. Indeed, up to a diagonal term 
the interaction Hj f j c can be written in the form: 

H JfJc ~J f _J + + J f+ J_ (5.10) 

with J denoting the total angular momentum. Since the core projected state is a super- 
position of components where the core function factor is a Wigner function of even values 
for K, the projected state characterizing the whole system is a superposition of compo- 
nents having a Wigner function with odd K as one of the factor states. The action of the 
raising or lowering operator on D J MK will transform it to D J MK+l or D J MK _ l with an even 
projection of J on the symmetry axis. The factor states associated to quasiparticles will be 
transformed from a K = 1 state either to a K = state or to a K = 2 state. Concluding, 
in the intrinsic frame the Coriolis term connects the K = 1 two qp states with the K — 
two qp states. In the case of odd particle-core system with one particle out of the core, 
the matrix element of the Coriolis interaction is different from zero only for the states with 
K = ±|. The interaction is attractive or repulsive depending on whether J + 1/2 is even 
or odd. Note that, by contrast, here a set of even number of nucleons is moving outside a 
phenomenological core and, moreover, the Coriolis interaction is effective in any iT-state. In 
the laboratory frame the spin-spin interaction affects mainly the 2qp states, as is shown in 
FigJHl It is worth noticing that the interaction is attractive in the 2qp band states having 
an angular momentum smaller than the crossing point and repulsive in other states of the 
S-band. Therefore, the spin-spin interaction causes an attenuation effect on the moment of 
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inertia in the S'-band, after crossing the g-band. If we switch off the spin-spin interaction, 
the slope of the curves from Fig. 4, representing the energy versus J (J + 1), is decreased 
after the two bands crossing point, which results in enlarging the moment of inertia of the 
S band. Due to this feature the bending in the moment of inertia plot is more pronounced 
than the bending taking place in the presence of the spin-spin interaction. 

The alignment of the particle angular momenta to the core angular momenta is shown 
in Fig. 7 b). The fact that the total K for the two neutrons is equal to unity prevents a full 

— * 

alignment. The larger J/ the larger the effect of the spin-spin interaction term. This infers 
that the backbending occurs when the broken pair is from a large spin single-particle state. 
For this reason the most favorable candidates for decoupling are particles from intruder 
orbitals like 213/2 for neutrons and hu/2 for protons. 

We stress on the fact that breaking a pair means to break the time reversal symmetry of 
the system, i.e to promote a paired particle to a state of different k which results in having 
a pair with a total K different from zero. Such an operation cannot be performed by the qQ 
interaction (see Eq.(3.15)) since the quasiparticle factor has only two quasiparticle K = 
terms. Therefore in the present formalism the only term responsible for pair breaking is the 
spin-spin interaction. 

Eqs. (5.7) and (5.8) can be used for calculating the gyromagnetic factor for an yrast state 
of angular momentum J. Indeed, from the expression of the magnetic moment 

a* = g c J c + 9fJf = gjJ, (5.11) 

one easily derive the following expression for gy. 

\ 4- ~\\ _ j ( T 4- ~\\ 

(5.12) 



93 = 9c + — ~ — 



J f (J f + l)-J e (J e + l) 
J(J + 1) 



where g c and gf denote the gyromagnetic factor of the core and fermionic system. For the 
core we consider 

9c ~ j-, (5.13) 

with Z c and A c denoting the charge and atomic number characterizing the core. Taking into 
account that the intruder state is 213/2 and this is occupied by neutrons, the corresponding 
gyromagnetic factor is: 

n 

9f = 13 ~ -°- 2943 ^' ( 5 - 14 ) 
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FIG. 7: The coupling scheme for g-band (a) and for 5-band (b) associated to a broken pair of 
particles which give rise to a Jj angular momentum and a general K = k\ + hi projection. The 
angular momentum of the particles is then coupled to that of the prolate core which is perpendicular 
to the symmetry axis (Z), resulting in the total angular momentum J. 

with \i N being the nuclear magneton and g™ the gyromagnetic factor for the neutron spin. 
The gyromagnetic factor plotted in Fig. 9 as a function of J, reflects the nature of the yrast 
states. Indeed, before the intersection of Oqp and 2qp bands, gj is independent of J and 
very close to the rotational value Z/A. When the intersection of the two bands takes place, 
gj has a big jump to a negative value, which confirms the 2qp character of the band which 
follows. This feature persists only for a few states and then the core contribution starts to 
be dominant which results in having a positive value for gj. The curve allure suggests a 
quadratic J dependence for gj. Except for 160 Yb, the curve has two branches, one constant 
corresponding to the g-band and one quadratically increasing with J, which corresponds 
to the S'-band. In the case of 160 Yb, after the crossing point, the curve decreases a little 
and then starts increasing. Note that because the magnetic moment carried by the core 
and fermions have opposite orientations and moreover the core contribution is an increasing 
function of J, there is a critical value of J where the total magnetic moment and therefore 
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FIG. 8: (Color online) The matrix elements of the spin-spin interaction are plotted as function of 
the total angular momentum. 

the gyromagnetic factor are vanishing. For the maximum angular momentum considered 
here, the gyromagnetic factor reaches half of the rotational value, i.e. about 0.2 fi^- 

Before closing this section we would like to spend few lines for comparing our approach 
with two other formalisms. At a superficial glance one may think that the formalism pre- 



sented here is similar to that of Ref. 
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45] . Therefore a fair comparison of the two proce- 



dures is necessary. First we mention that the mean field determining the single-particle basis 
is deformed while that used in Ref. 
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45j is spherical. Therefore, an angular momentum 



projection operation for the deformed BCS states is necessary in our case. In the quoted 
reference, the coherent state is defined in terms of deformed quadrupole bosons which makes 
difficult the evaluation of the overlap matrix elements. Indeed, the matrix elements between 
rotated intrinsic states are approximated and then parametrized. By contrast, the coher- 
ent state used here corresponds to spherical quadrupole bosons and moreover the matrix 
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FIG. 9: (Color online) The gyromagnetic factor given in units of nuclear magneton, calculated for 
yrast states, is represented as function of angular momentum. 



elements are analytically calculated. Here the pairing interaction is effective for nucleons 
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only the pairs from the core 



moving outside the phenomenological core, while in Ref. 
and those from outside the core interact with each other. Moreover, the pairing interaction 
is treated in the particle representation which in fact leads to a parametrization of the cor- 
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responding matrix elements. Consequently, the Hamiltonian is diagonalized in a basis with 
a fixed number of particles. Since a BCS formalism for the deformed single-particle states is 
used, in the present case, the number of particles is conserved only in average. The numbers 



of the parameters employed by the two formalisms are also different: seven for Ref. |31l. 145 1 
and five in the present approach. One may conclude that although the two approaches have 
some common features, they are essentially different. A similar single-particle basis is used 



in Ref. 



46| but in a different context. Indeed, using a particle-rotor formalism the depen- 



dence of the interaction of the lowest two bands on the degree of filling the shell is studied 
pointing out an oscillating behavior. 



VI. CONCLUSIONS 



In the previous sections a semi-phenomenological formalism for the description of the 
backbending phenomenon was proposed. A model Hamiltonian associated to a set of in- 
teracting particles moving in a deformed mean field coupled to a phenomenological core 
described in terms of quadrupole boson operators is treated in a product space of angular 
momentum projected states. The pairing interaction of neutrons moving in a deformed 
mean field is treated by the BCS formalism. The model states for the ground state band 
are obtained by angular momentum projection of the deformed product state \BCS)di> c , 
while the S'-band states are projected out from the intrinsic K = 1 two quasiparticle states 
J + aJ fc aj_ fc |£>C S) dtpc- The sub-state \jk) is chosen such that the corresponding quasipar- 
ticle energy is minimum. Projected states of g and S bands are not mutually orthogonal. 
Diagonalizing the overlap matrix, one defines an orthogonal basis for treating the model 
Hamiltonian. The parameters involved were fixed by a fitting procedure described in the 
previous section. The lowest Hamiltonian eigenvalues in the orthogonal basis defines the 
yrast band. The first energy levels originate from the projected states of Oqp state mentioned 
above while starting from a critical angular momentum they are mainly of 2qp nature. 

The experimental backbending shape of the moment of inertia versus angular frequency 
squared is fairly well reproduced by our results. The data considered in our calculations 
refer to energy levels with angular momentum up to 26-28. A measure of the agreement 
quality is the r.m.s. value characterizing the deviations of the calculated energies from the 
corresponding experimental data. This quantity is about 30keV or less. 
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A detailed analysis of the effect coming from the Hj f j c term concerning the neutron 
pair breaking as well as the alignment of individual angular momentum to the core angular 
momentum, is presented. The pair breaking takes place for J = 10, 12 and the maximum 
alignment is settled for J larger than 20. Since the de-paired neutrons carry a projection 
K — 1, a full alignment is not possible. As shown in Fig. 5 for J > 20 the angular momentum 
defect AJ reaches a plateau with AJ = 2. It is also observed an abrupt change of the 
gyromagnetic factor from positive to negative values in the band crossing region, which 
is consistent to the change in structure of the yrast band, going from rotational to 2qp 
character. 

We may conclude that the present formalism is able to account quantitatively for the 
main features of the backbending phenomenon in the rare earth nuclei. 

Although effects like pair breaking, bands crossing and angular momentum alignment 
have been explained, within some cranking formalisms, by several authors, the present 
approach provides a consistent description for the mentioned effects pointing out several 
new features which will be enumerated below. 

• Although we use a spherical projected particle-core basis the two components, particles 
and core, are deformed. Working with states of good angular momenta is an advantage 
over the cranking methods where the states in the crossing region exhibit a large 
dispersion for angular momentum. 

• The effects of the qQ and Jf ■ J c terms on the backbending are discussed and specific 
contributions are identified. 

• It is shown that the spin-spin interaction simulates in the laboratory frame the Coriolis 
force which is active in the intrinsic frame. 

• By the bands crossing the core angular momentum, J c , varies only by a few units, 
2 — Ah. The minimal variation is of 2h and is recorded for 160 Yb. Due to this feature 
one expects a nonvanishing E2 transition between the states lying in the vicinity of 
the crossing point. 

• We depicted one case, 160 Yb, where after the bands crossing point, the angular mo- 
menta carried by fermions and the core are close to each other. 
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• The curves associated to J c and J/, in Fig. 5, cross each other at a total angular 
momentum equal to 22 for 156 Dy, 158 Er, 24 for 160 Yb, 160 Er and 26 for iHi66 Hf _ 

• The maximum alignment is reached in the region of J c > J/, where the defect A J is 
equal to 2 — 3h. This is a reflection of the K — 1 nature for the S-band. 

Due to the above mentioned aspects one may say that although the present paper ad- 
dresses a relatively old subject the proposed formalism unveils new features of the back- 
bending phenomenon. 

Before closing, we present few perspectives of the present approach. Indeed, the results 
encourage us to extend the restricted model space by adding to the particle factor the proton 
state /in/2 which is suspected to be responsible for the second backbending. Another possible 
extension refers to the collective factor states, by adding to the coherent state considered 
here the model states for beta and gamma bands used by the coherent state model. It is 
well known that the energy spectra of these bands comprise more irregularities than the 
ground state band. It is an open question whether such anomalies could be also interpreted 
as the interaction with other bands of different nature. In this way we could describe a 
multibackbending phenomenon showing up in the yrast and non-yrast bands. 
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VII. APPENDIX A 

The analytical expressions for the terms Ai(M, J) and B n (M,J) involved in Eq.(2.10) 



are: 



M 




7Y 



(A.l) 



m ■ n 



M + 1 




(A.2) 



A(M, J) 



(A.3) 



UL (M-J + 2k+l) 
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The integral Ib(J,m) can be analytically determined, and the result for the present case 
J =even and m =even is: 

f 2 nfc^ 2 -(2A:) 2 ] m > Jf 
I b (J, m) = I Yll =0 [m 2 - (2k + l) 2 ] ' (A.4) 
[ 0, m < J. 

VIII. APPENDIX B 

The state obtained by applying a pair of time reversed quasiparticle operators applied 
on a BCS function can be written as a linear combination of states with a definite angular 
momentum J, with J running from to 24, and a projection on z-axis K — 0: 

a) k a)_ k \BCS) d = Y. C f\ J ^)- 
j 

Acting with an angular momentum projection operator on this state is equivalent to selection 
of only one J component from the above linear combination, and rotate the projection to 
the value M: 

P J J f AA-k\BCS) d = C*\J f , M f ). (B.2) 

The angular momentum states | J, M) are orthogonal and normalized to unity. Then the 
norm of this projected function is the reciprocal of the amplitude C^: 

{C^) 2 = d (BCS\a_ k a k P^a\al k \BCS) d . (B.3) 
Now if we apply on the state (B.l) a raising angular momentum operator, 

J + a) k a)_ k \BCS) d = £cJV^>l)l4 1), (B-4) 

j 

and then project an angular momentum from the resulting function, one obtains exactly the 
K — 1 2qp projected fermionic function: 

P'M f iJ + ^U\BCS} d = C>j k f y/j f (J f + l)\J f ,M f ), (B.5) 

whose normalization factor is readily obtained: 

(At,) " = Wf + 1) d (BCS\a- k a k P^aW_ k \BCS) d . (B.6) 
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Making use of Eqs. (B.5) and (B.6), one can also derive the expressions for other m.e. 
which we need in our calculations. For example the overlap of Oqp and 2qp projected states 
can be expressed as: 



d (BCS\P^J^U\BCS) d =[N^N^ l \ ■ (B.7) 

In order to calculate the quantities d {BC S\aj_ k aj k PQQ a^ k ,a^_ k ,\BC S) d first are calcu- 
lated the corresponding matrix elements with particle operators Cj k (cj k ) instead of quasi- 
particle ones otj k {otjk)i for which we use the same method as in the case of the norm of 
the projected BCS state. Having these determined, the desired matrix elements are easily 
obtained by applying a canonical transformation to the particle operators: 

d (BCS\a j - k a jk P^aLa]_ k ,\BCS) d = \ d (BCS\ Cj - k c jk P^ cLc]_ k ,\BCS) d 



jk' I / " TT2 TV1 

u jk u jk' 



■^^(-) k+k ' d (BCS\P^\BCS) d - j^{-r k d {BCS\P J ic] k ,c)_ k ,\BCS) d 



^—{-y- k ' d (BCS\ Cj - k c jk P^\BCS) d , (B.8) 



U 3k' U jk 

where U and V are the occupation parameters from the BCS equations. 
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